Ginzburg-Landau equation bound to the metal-dielectric interface and transverse 
nonlinear optics with amplified plasmon polaritons 
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Using a multiple- scale asymptotic approach, we have derived the complex cubic Ginzburg- 
Landau equation for amplified and nonlinearly saturated surface plasmon polaritons propagating 
and diffracting along a metal-dielectric interface. An important feature of our method is that it 
explicitly accounts for nonlinear terms in the boundary conditions, which are critical for a cor- 
rect description of nonlinear surface waves. Using our model we have analyzed filament at ion and 
discussed bright and dark spatially localized structures of plasmons. 



I. INTRODUCTION 

Surface plasmon polaritons (SPPs) are half-photon 
half-electron surface waves. Thanks to their dual na- 
ture, SPPs can be focused tighter than pure Hght, which 
is an important property for potential applications in op- 
tical processing of information. In the absence of lateral 
boundaries propagating SPPs are expected to diffract 
in the interface plane. One way to control diffraction 
is to structure the surface and make plasmonic waveg- 
uides, see, e.g., [1]. Alternatively, one can use nonlin- 
earity for the creation of spatial SPP solitons, which are 
non-diffracting self-localized surface waves (2l,[3(. Other 
transverse effects with nonlinear plasmons, such as self- 
focusing and filamentation [4] can be important for fre- 
quency conversion, switching and routing experiments 
with photonics chips. Note, that the interplay between 
transverse and nonlinear effects have attracted significant 
attention outside the plasmonics and nanophotonics con- 
texts, see, e.g., [5|, |6|] for the historic accounts. 

Nonlinear functionality of SPPs 0-[2| can be signifi- 
cantly hampered by ohmic losses resulting in short prop- 
agation distances. One of the possible solutions is to 
amplify SPPs by doping and pumping the dielectric, so 
that the losses are either partially or fully compensated 
p!Ql-[l2j. The linear dispersion of the amplified SPPs has 
been studied by several groups, see, e.g., [13, 14]. The 
linear results have been recently generalized to the more 
realistic case, when linear gain is nonlinearly saturated 
above the stimulated emission threshold [15]. 

Analytical or semi-analytical approaches to describe 
nonlinear effects with SPPs are very important, since 
the first principle numerical modelling of nonlinear and 
multidimensional cases is still computationally demand- 
ing. Recently the nonlinear Shrodinger equation (NLS) 
has been introduced for the plasmons in a slot waveguide 
formed by two planar metal dielectric interfaces j2| and at 
a single interface [3[ . The averaging method implemented 
in [2|,l3[ has been borrowed from the theory of the dielec- 
tric waveguides [l6|. In this approach, one starts from 
the known solution for the linear SPPs: F(x)e*^^, where 



X is the coordinate perpendicular to the interface, z is the 
propagation direction and j3 is the propagation constant. 
Then introducing a slowly varying amplitude A{z., y) and 
assuming small nonlinearity, the Maxwell equations are 
averaged in x and the NLS equation for A is derived @,|3[. 

The above approach has some drawbacks. First, it 
is sufficiently well justified only for the quasi-transverse 
fields, that approximately satisfy the wave equation [l6|. 
Another problem is that it does not treat the boundary 
conditions rigorously. In particular, continuity of the nor- 
mal to the interface component of the displacement Dx 
is guaranteed only in the linear approximation. If the in- 
tensity of the guided light peaks away from the interfaces 
and is small in the proximity of the latter (like it typically 
happens in dielectric waveguides operating on the princi- 
ple of total internal reflection), then the nonlinear terms 
in the boundary conditions can be disregarded. However, 
the SPP intensity peaks exactly at the interface and non- 
linear contribution to the boundary conditions is critical 
[l5|, [13]. This is also true for other types of nonlinear 
surface waves, see, e.g., [18|-[2Q|- 

Thus, it is important to develop a rigorous procedure of 
deriving the nonlinear evolution equation for surface plas- 
mons, such that the continuity of is enforced together 
with its nonlinear part. Below we develop a multiple- 
scale asymptotic approach for the amplified SPPs, which 
treats the nonlinear boundary conditions rigorously and 
reveals the differences with the results obtained by the 
averaging technique. The gain and complex nonlinear- 
ity we use are derived using the two-level model. Our 
procedure leads to the complex Ginzburg-Landau equa- 
tion that accounts for diffraction of SPPs in the interface 
plane. The nonlinearity enhancement factor derived by 
us is intrinsically complex, whilst the averaging approach 
gives a real one. The difference between the predictions 
of two approaches increases in the short wavelength limit, 
where SPPs are maximally localized at the interface, 
and thus the nonlinear part of the boundary condition 
is more important. Using our theory we derive a crite- 
rion for SPP filamentation and discuss bright and dark 
spatially localized SPPs. Strictly speaking both soliton 
families are unstable, though the bright solitons demon- 
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strate propagation distances sufficient for their practical 
observation. 



II. MODEL 

We assume that the interface between metal and dielec- 
tric is at X = and y are the in-plane coordinates. The 
evolution of SPPs obeys the time independent Maxwell 
equations 

^ly^y ~ ^yy^x - d^z^x + d'^x^z = e^x, (1) 

9yz^z - dl^Ey - dl^Ey + dlyEx = eEy, (2) 
^Iz^x - dlx^z - dyyEz + d'^yEy = eEz. (3) 

The coordinates are dimensionless and normalized to the 
inverse wavenumber k = 27r/A^ac, where Xyac is the vac- 
uum wavelength. The permittivity on the dielectric side 
of the interface {x > 0) is 

e = ed^Xmf^\Eyf^\Ezf), (4) 
= 4 + ^4'. X = X ^ ^X'- (5) 

The permittivity on the metal side {x < 0) is 

If SPPs are amplified by means of active inclusions in 
the dielectric, then Cd and x functions of the gain 
coefficient a. The propagation constant /3 for Hnear plas- 
mons is 

P V^d^m/{^d + ^m)' (7) 

where /3 becomes real at the threshold: a = ao [l3|: 

P{ao) = Po. ImPo = 0. (8) 

The linear and nonlinear permittivities for the dielectric 
at a = ao are ed{ao) = Cdo and x(<^o) = Xo- 

The active inclusions are approximated by the 2-level 
atom susceptibility. For light intensities much smaller 
than the transition saturation intensity Ig we find j^, I^H 

i — S 

Cd = eb- c^Y^rp' 

X = ^(i^- (10) 

where is the real dielectric constant of the background 
material hosting the two- level atoms. S = {uj — uJa)T2 
is the dimensionless detuning from the atomic resonance 
frequency, uJa = 27rc/Aa, normalized to the transition 
linewidth, T^^. a is the dimensionless gain coefficient 
at the line center. The electric field is normalized to 
V^g, which implies that the nonlinear susceptibility x 
is dimensionless, see [15| for more details. Possible de- 
pendence of the atomic life times from the distance to 



the interface, see, e.g., [22*1, are specific to a choice of a 
pumping technique and are disregarded in what follows. 
The threshold gain ao works out as 




Lossless propagation of SPPs is impossible above the crit- 
ical value of S = Sum as determined by the condition 
that the square root in Eq. (pT|) becomes zero. At this 
point, the two solutions for ao degenerate, see Fig. [H 
Slim should not be confused with the plasmon resonance 
frequency, Sspp^ which corresponds to the zero of the de- 
nominator of Po' The existence boundaries of the SPPs 
dit a = ao are determined by either or both of dum and 
Sspp^ see Fig. [H The upper branch solution (dashed lines 
in Fig. [1]) corresponds to high gain coefficients imply- 
ing refractive indices of order 10 or larger. In our sub- 
sequent numerical examples we focus on relatively small 
|(5|'s, thereby selecting the lower branch of ao (minus sign 
in front of the square root in Eq. (fTT]) ). This branch cor- 
responds to relatively small changes of the background 
refractive index as achievable for small densities of ac- 
tive atoms [14]. Note, that Cm is frequency dependent, 
i.e., em(^) = em{S/T2 +cja)- Hence ao is a function of 
both 6 and uJa- We choose silver as a metal in all our 
calculations. 



III. FIRST PRINCIPLE DERIVATION OF THE 
GINZBURG-LANDAU EQUATION FOR SPPS 

The perturbation theory developed below assumes rel- 
atively small deviations of the gain coefficient from its 
threshold value ao, i.e. 

a - ao = ao^, \g\ <C 1. (12) 

The anzats for the field components in the dielectric and 



metal is 
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j = d,m 


where A^^] ^ 




x,j 







and 4^] ~ |5|3/2. All A's are the functions of 
y and x. However, their dependencies on z and y are 
assumed slow relative to the fast oscillations of e*^°^: 



dy - \g\"^ ~ \9\. (14) 

Though the y component deviates from zero, when the 
field has finite size along it is expected to remain rel- 
atively small, ~ |^|. 
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FIG. 1: (Color online) (a) Threshold gain ao vs detuning b for 
two different atomic resonances: \a — 400nm (line 1 (blue)) 
and \a — 700nm (line 2 (red)), (b) /3o vs 5. Parameters and 
notations as in a). The full/dashed line corresponds to the 
minus/plus in Eq. (pT]) . 



(18) 



Any nonlinear and transverse, i.e., ^-dependent, effects 
are disregarded in Eq. (p!6]) . 

The SPP solution of Eq. ([I6j) is well known: 



.(0) 
j(0) 

^(0) 



Id 

Qm 



(19) 



Eqs. (p!9]) satisfy continuity of the normal component 
of the displacement and tangential components of the 
field: CdA^^^^ = CmA^^^m and = at x = 0. The 

former condition implies edoQm = —^mQd giving (after 
some algebra) the expression for /3o, see Eq. ([7|). 

In the order |^|^, we find the linear equations for the y 
component of the field 

i(0) 



which are readily solved 
^(0) _ 

^(0) = 



8^ 



0, 



5(y,z)e-^^^ 



(20) 
(21) 



To determine the unknown functions A{y^z) {\A\ ~ 
l^l"*-/^) and B{y^z) {\B\ ~ |^|) we need to proceed to 
the higher orders of our perturbation series. 



.|3/2 



order and Ginzburg-Landau equation 



The dielectric susceptibility is trivially expanded into 
the ^f-series: 



= ^dO + ^dlj ^dl = gOL^da^d- 



(15) 



The only value of the nonlinear coefficient Xi we require 
below is the one taken exactly at the threshold x(<^o) = 
Xo- 



A. \g\^^'^ and l^^l^ orders 

By substituting Eqs. ([T3|) into the Maxwell equations, 
we find in the order \g\^^'^ 



A 



i(0) 

(0) 



= 0, j = d,m. 



Here 



9] 




(16) 



(17) 



Proceeding to the order we find an inhomoge- 

neous system of differential equations for corrections to 
the standard SPP solutions. The correction equations on 
the metal side are 



Mrr 



'^x,m 













(22) 



where 



^l±^d,A-q^dyB-'-^dlA, 

Qm Qm 



-2i/3od.A - dl^A. 



Kx 

A solution of Eqs. (|^ consists of a particular solution 
of the inhomogeneous problem plus a general solution of 
the corresponding homogeneous system {K^^z = 0): 

1 



— c — e^'" 

Qm 

ZQm 



(23) 
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where c is a constant to be determined from the boundary 
conditions. 

The righthand sides of the corresponding equations 
in the dielectric are more cumbersome due to nonlinear 
terms: 



(1) 

(1) 
z,d 



(24) 
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where 



Id 



ij3o 

qd 



ediA- 



+qadyB , 
L, = -2ipod.A - ediA - d^yA 
N. = '^(j% + l)xo\AU 

qd Vkdl / 



(erfo^d + 2/?2Re(Zrf) {\qd\^ + Pi] 



Solutions of Eqs. 



^doqdmr 

are given by 



■Xo\A\'A. 



1 



ipoN,{2Reqd^qdy 



^Qd 



(25) 



d L 



Nx 



4Reqd{Reqd ^ qd) _ 



-2ReqdX — qdX 



-qdX 



2qd 4Reqd(Reqd ^ qd) 



The arbitrary constant terms in Eqs. (|25]) have been 
omitted, as this does not lead to any loss of generality. 

Combining Eqs. ([I9|) , (|2l]) , (|23|) , ((S^ (p^ and 

substituting the calculated fields into the boundary con- 
ditions 



-^x,d 



[ed^x{\Ex,d\^^\Ey,d\^^\Ez,d 

— ^m-E^x,mi (^^) 

we find that the latter are satisfied in the order \g\^^'^ 
only providing that 



^qd + Reqd)Reqd 



(27) 



and the amphtude A solves the complex Ginzburg- 
Landau equation 



2il3od,A + dlyA + fA + 7|A|2^ = 0, 



(28) 



where 



7 /ixo, h 



(qo + em)^' 

P'o qd{\qd\'^Pi) 



^do (qd^^eqd) \qdf 



All the terms containing B{y^ z) cancel out leaving this 
function undetermined until the higher order corrections 
are accounted for. Thus taking the plasmonic field as in 
Eqs. (p!3|) with the amplitude A obeying Eq. (|28|) we are 
guaranteed that the nonlinear boundary conditions are 
satisfied upto and including the |^p/^-terms. 

The first and second terms in Eq. (|28|) describe the 
propagation and diffraction of SPPs. Re/ accounts for 
the shift of the propagation constant away from /3o, 
when gain deviates from the threshold. Imf accounts 
for the gain excess {a > ao, Imf < 0) or shortage 
{a < Qfo, Imf > 0). The nonlinear term provides an 
additional shift of the propagation constant (Re7|Ap) 
and of the nonlinear loss (Im7|Ap) that counterbalances 
the linear gain. Note, that the transformation back to 
physical units results in the appearance of a factor in 
the 3rd and 4th terms of Eq. 



C. Comparison with the averaging approach 

Parameter h in the expression for 7 is the nonlinearity 
enhancement factor, h accounts for the difference be- 
tween the nonlinear responses of SPPs and free waves 
propagating far from the interface, h is complex, there- 
fore, even when atomic nonlinearity is purely dispersive 
or purely absorptive, the effective SPP nonlinearity is a 
mixture of both types. This contrasts results using the 
averaging approach [l, [sl, [16] , where h is real. 

The averaging approach yields a well known expression 
for the effective nonlinearity of guided modes. Following 
this method, one should replace h in Eq. (l28]) with h 



Jo 



J —OD 



(29) 



Here F is the plasmonic field given by Eq. (fT9|) with 
A=l. 

We fix the detuning and plot h and h as functions of 
the resonance wavelength, in Fig. [2l On the short 
wavelength side the plots in Fig. [2] are limited by the 
zero of the denominator of Pq. Through this it is seen 
that the two approaches give qualitatively similar de- 
pendencies in the long wavelength limit, whilst in the 
short wavelength limit, our calculations predict a signifi- 
cantly higher nonlinearity enhancement. Physically, one 
can identify two factors determining changes in h with 
the resonance wavelength. Tendency for Reh and h to 
decrease with decreasing is linked to the fact that 
SPP intensity on the metal side increases relative to the 
intensity on the dielectric side, as wavelength decreases. 
Since the metal is linear in our model, it should lead to a 
drop in the nonlinearity enhancement coefficient. How- 
ever the smaller wavelengths reaching the SPP resonance 
imply that the SPP field profile is getting squeezed closer 
to the interface on both sides and therefore the nonlinear 
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FIG. 2: (Color online) (a) Nonlinearity enhancement coeffi- 
cients Re/i (full line) and h (dashed line) calculated using two 
different approaches vs Xa. (b) Imh. Varying S inside the 
transition linewidth leads only to small difference in h. The 
graphs shown correspond to 6 = —0.5. The short wavelength 
boundary of the both plots corresponds to the point where /3o 
becomes imaginary. 



part of the boundary conditions becomes more impor- 
tant. This makes the difference between Keh and /i, and 
the deviation of Imh from zero pronounced in the short 
wavelength limit, see Fig. [2l 



IV. FILAMENTATION AND TRANSVERSELY 
LOCALIZED SPPS 

A. Filamentation of SPPs 



The plane wave solution of Eq. 

z 



IS 



Aq = p X exp 
2 Im/ 



z--{Ref - p^Rej) 



(30) 



Im7 



> 0. 



— Im / > implies the gain above threshold, i.e., a > a^. 
Im7 > implies absorptive nonlinearity compensating 
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FIG. 3: (Color online) Growth rate k of the filamentation 
instability expressed in physical units as a function of mo- 
mentum p. Xa — 594 nm, S — —0.3, ao — 0.0063 and 
g = 0.4,0.5,0.6 for blue (line 1), green (line 2) and red (line 
3) curves respectively. 



for the excess gain. Together these conditions lead to 
the existence of the SPP solution with the constant sta- 
tionary amplitude, p. The expression for p, however, 
diverges and our approach breaks down if the nonlinear 
absorption becomes zero, i.e. Im7 = 0. This critical case 
requires one to account for quintic nonlinear terms in the 
perturbation expansion, which goes beyond our present 
objectives. 

The solution (|3Q]) can be unstable with respect to a 
pattern forming filamentation instability, as known for 
the generic Ginzburg-Landau equation [23|. In order to 
study the stability problem we perturb Aq with small 
amplitude waves carrying transverse momentum p: 



A = (1 + g+e 



Kz-\-ipy 



q_e 



K z—ipy 



)Ao. (31) 



Inserting Eq.(j3Tj) into Eq. ([28|) and linearizing for small 
\q±\ we find two solutions for The unstable one is 
given by: 



2pon = Im/ + ^ImV - pHp^ - 2p^„J. (32) 

The filamentation instability sets in providing Kehc > 0. 
In Fig. [3l it is seen that the Rehz vs p plot has the typical 
two peak shape. The maximal instability growth rate is 
achieved for 



p = ±pr, 



Pr, 



(33) 



The characteristic filament size in physical units is w ^ 
^vac/Pmax' w diS di fuuctiou of is plotted in Fig. |4l 

The instability domain in the (J, a)-plane is shown in 
Fig. [5l Filamentation is present for the self-focusing 
effective nonlinearity, i.e., Re7 > 0. If the nonlinear- 
ity enhancement factor is real, then Re7 = simply 
implies Imx = 0, which is achieved at the line center 
6 = 0. In this case, the nonlinearity changes from focus- 
ing (filamentation) to defocusing (no filamentation) at 
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FIG. 4: (Color online) The characteristic filament size w 
scaled back into physical units vs Xa for 6 = —0.5, g = 
0.1,0.3,0.5 for blue (line 3), green (line 2) and red (line 1) 
curves respectively. 




FIG. 5: Lossless SPPs exist above the full line corresponding 
to a = ao- SPPs are unstable with respect to filamentation 
on the left from the dashed vertical line, corresponding to 
Re7 = 0. Aa 400 nm. 



the atomic resonance, i.e. exactly as in the bulk mate- 
rial. However, the fact that Imh ^ leads to the offset 
of the instability boundary away from ^ = 0, see Fig. [5l 
Gain values corresponding to approximately 50% above 
threshold imply the development of a filamentary pattern 
over distances of 1 — 3mm. 



B. Bright and dark localized SPPs 

The cubic Ginzburg-Landau equation is known to have 
a wide variety of localised solutions, which can be rele- 
vant in the SPP context under different circumstances. 
Detailed classification of these solutions can be found in, 
e.g., Refs. j23|, [IJ. Here we briefly introduce the most 
ubiquitous of those, which are bright Pereira-Stenflo [2^ 
and dark Nozaki-Bekki i26il localised solutions. 



FIG. 6: (Color online) Propagation of a bright localized 
solution and instability of its background: \a — 594 nm, 
b = -0.3, OL = 0.0183, ao = 0.0063. 
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FIG. 7: (Color online) Destabilization of the dark soliton due 
to core instability: \a — 594 nm, b — —0.3, ot — 0.0183, 
ao = 0.0063. 



Both bright and dark localised solutions exist under 
the same conditions: Im/ < (positive gain) and Im7 > 
(nonlinear absorption) . The bright solution is given by 



and the dark one by 



(34) 



. tanh(sy) , . . 
[cosh{sy)\ 

Explicit expressions for the parameters entering Eqs. 
(|34|) . (j35j) are given in the Appendix. In the limit 
\y\ ^ oo the dark solution tends towards the plane wave 
solution (|3Ql) . 

The bright solution is unstable because its zero back- 
ground is unstable above threshold. This instability is 
relatively slow to develop and practical observation of 
bright solitons over distances of 100s of microns is still 
possible, see numerical modelling results in Fig. [6l The 
dark solution is known to be unstable with respect to the 
core instability through most of its existence domain, see, 
e.g., which is complemented by the filamentation of 
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background, provided the effective Kerr nonlinearity is 
self- focusing. Fig. [3 shows an example of the core insta- 
bility, where one can see that for the chosen parameters 
it develops over the shorter distance, if compared to the 
instability of the zero background of the bright solution 
in Fig. [6l 



SUMMARY 



anisms leading to stabilization of these structures is an 
important topic for future research. 
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We have considered nonlinear propagation of the am- 
plified and diffracting surface plasmon polaritons above 
the threshold beyond which the plasmon propagation 
constant becomes real. Starting from the first principle 
Maxwell equations, we have developed a technique al- 
lowing derivation of the complex cubic Ginzburg-Landau 
equation for the slowly varying plasmon amplitude. The 
nonlinear plasmon solutions found by our method is guar- 
anteed to satisfy nonlinear boundary conditions at the in- 
terface to the required accuracy. This distinguishes our 
approach and results from the recently proposed deriva- 
tion of the nonlinear Schrodinger equation for surface 
plasmons, which satisfy boundary conditions only in the 
linear approximation [l, Q . We have found that the non- 
linearity enhancement factor is always complex and hence 
mixes real and imaginary parts of the intrinsic nonlin- 
earity of a dielectric. This mixing changes conditions 
required for the filamentation of plasmons and the ex- 
istence of the dark (Nozaki-Bekki) and bright (Pereira- 
Stenflo) spatially localized waves, relative to the respec- 
tive conditions in bulk medium. Though both of the 
localized solutions are unstable with respect to growth of 
small perturbations, the bright ones demonstrate quasi- 
stable propagation over the distances of lOO's of microns 
and are likely to be practically observable. Finding mech- 



Appendix 

The parameters entering Eq. (|34|) are expressed in 
terms of the parameters for the Ginzburg-Landau equa- 
tion Eq. dSHl) as 




1 - 1 

2/3o 4/3oa 
Parameters entering Eq. ([35]) are 



2Im7 



/3Re7y 
V2Im7 J 



3b 

V = — — Ref — — Imf. 
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